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Abstract 

Modern physics is largely devoted to study conservation laws, such as charge, energy, linear 
momentum or angular momentum, because they give us information about the symmetries of 
our universe. Here, we propose to add the relationship between electromagnetic duality and he- 
licity to the toolkit. Generalized electromagnetic duality symmetry, broken in the microscopic 
Maxwell's equations by the empirical absence of magnetic charges, can be restored for the macro- 
scopic Maxwell's equations. The restoration of this symmetry is shown to be independent of the 
geometry of the problem. These results provide a simple and powerful tool for the study of light- 
matter interactions within the framework of symmetries and conservation laws. We apply such 
framework to the experimental investigation of helicity transformations in cylindrical nanoaper- 
tures, and we find that the transformation is significantly enhanced by the coupling to surface 
modes, where electromagnetic duality is strongly broken. 



I. INTRODUCTION 



Symmetries, both continuous and discrete, are a powerful tool for studying Nature. Ac- 
cording to Noether's celebrated theorem PQ, any continuous symmetry of a non-dissipative 
system gives rise to a conserved quantity in the dynamic equations. In modern algebraic 
terms we say that when a system is invariant under the continuous transformation generated 
by a given operator, the observable represented by that operator is a conserved quantity. For 
example, rotational and translational invariance are associated with conservation of angular 
momentum and linear momentum because, as transformations, rotations are generated by 
the components of angular momentum and translations are generated by the components of 
linear momentum. 

After Einstein founded his theory of space and time on relativistic invariance, the study 
of symmetries has never left center stage in modern physics. Wigner found that fundamental 
particles are defined as objects that are invariant under certain transformations: The same 
principle is currently used for both standard model and beyond standard model physics. 
Many of the symmetry transformations used in physics have a geometrical nature: trans- 
lations, rotations, Lorentz boosts. Some others have a more abstract character, separated 
from geometry, like the isospin symmetry which groups together protons and neutrons into 
nucleons or the color symmetry of hadrons, which governs the strong force interactions. In 
light matter interactions, symmetry considerations allow for instance to establish rules for 
exciting atomic electrons with light. Symmetry reasons are also behind the few exact solu- 
tions of Maxwell's equations in inhomogeneous media, like planar multilayers and spheres. 
Other geometries, like finite cylinders, are also best studied by exploiting their symmetries, 
as will be shown in the scattering experiments off cylindrical nanoapertures that we present. 

In this paper we will study a non-geometrical symmetry in electromagnetism in detail: 
electromagnetic duality. Electromagnetic duality is a transformation where the roles of elec- 
tric and magnetic fields are mixed. Mathematically, the generalized duality transformation 
of electromagnetic fields is expressed as: 
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In the absence of charges and currents, ([T]) is a symmetry of Maxwell's equations: If the 
electromagnetic field (E, H) is a solution of the free space Maxwell equations, then the 
field (E e , H#) is also a solution. In 1965, Calkin [2] showed that helicity was the conserved 
quantity related to such symmetry. 

Helicity is defined [31 chap. 8.4.1] as the projection of the total angular momentum J onto 
the linear momentum direction P/|P|, i.e. A = J • P/|P|. In the case of photons jH chap. 
2.5], helicity takes the values ±1. It is possible to intuitively understand the meaning of 
helicity when considering the wave function of the particle in the momentum representation, 
that superposition of plane waves. In this representation, helicity is related to the 

handedness of the polarization of each and every plane wave. Only when all the plane waves 
have the same handedness is the helicity of the particle well defined. Note that polarization 
in momentum space and polarization in real space are not the same concept. What Calkin 
showed is that, as an operator, helicity generates generalized duality transformations in the 
same way that linear momentum generates translations and angular momentum generates 
rotations. Since that seminal work, the role of helicity as the generator of generalized duality 
symmetry transformations for the free space Maxwell's equations has been reported several 
times [SH?]. 

In 1968, Zwanzinger [5] extended this free space invariance and conservation law to a 
material quantum field theory with both electric and magnetic charges. In material systems, 
the symmetry studied by Zwanzinger seems to be broken by the lack of experimental proof 
of the existence of magnetic charges. The experimental efforts to find magnetic monopoles 
are a very active field of research [9HTT] because, as Dirac demonstrated [12], the mere 
existence of a single magnetic monopole in the universe would explain the quantization 
of electromagnetic charge. In Zwanziger's work, transformation Q is complemented with 
a corresponding mixing of electric and magnetic charges and currents. In the absence of 
magnetic charges, the microscopic Maxwell's equations are no longer invariant under the 
generalized duality transformations of ([!]). This is the current status of duality symmetry 
in material systems. 

In this article we show that the generalized electromagnetic duality symmetry, broken for 
the microscopic Maxwell's equations by the absence of magnetic charges, can be restored 
for the macroscopic Maxwell's equations for material systems characterized by electric per- 
mittivities and magnetic permeabilities. The restoration condition for a system composed 
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of different isotropic and homogeneous domains depends only on the materials and is inde- 
pendent of the shapes of the domains. When the system is dual, the helicity of the light 
interacting with it is preserved. 

The geometry independent extension to material systems presented in this article turns 
the relationship between helicity and duality into a simple and powerful tool for the practical 
study of light-matter interactions using symmetries and conserved quantities. The practical 
applicability of our ideas is enabled by a fact of crucial importance: Measurement and 
preparation of light beams with well defined helicity can be done with very simple optical 
elements. Armed with this tool, we experimentally investigate helicity transformations in 
focused light fields that interact with cylindrical nanoapertures in a gold film over a glass 
substrate. The study of the symmetries of the system allows us to identify the exact reason 
for the relatively large helicity conversion that we find in the nanoapertures: The coupling 
of light to surface modes, which, being equal weight superpositions of modes of opposite 
helicity, strongly break duality symmetry. This result shows the ability of the framework to 
make both qualitative and quantitative predictions. 



II. HELICITY AS THE GENERATOR OF GENERALIZED DUALITY TRANS- 
FORMATIONS IN FREE SPACE 

In our derivations, we will use a harmonic decomposition of the fields and assume a 
exp(— iut) dependency with the angular frequency u. Additionally, we will work in the 
representation of space dependent vectorial fields, also known as the real representation. 
This setting is different from those in [2] and [5], and, although the final result is not new, 
the derivation in this section sets the stage for the study of the piecewise homogeneous and 
isotropic case. 

The expression of the helicity operator for monochromatic fields in the real representation 
can be obtained directly from the definition of helicity: 

J^P _ (S + L)-P _ST_Vx 

|P| ~ |P| ~ |P| ~ k ' { ) 

where S and L are, respectively, the spin and orbital angular momentum operators, the 
third equality follows from the orthogonality of L = r x P and P, and the last one is valid 
in the real representation because S ■ P = Vx [131 expr. XIII. 93] and |P| is equal to the 
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wavenumber k for monochromatic fields. 

Related to the different settings mentioned above, a clarification regarding different def- 
initions of helicity is in order before we start. In [2J, helicity appears as an operator in the 
Fock space representation and in [8] as an integral involving the electric and magnetic fields 
and potentials operators. Both of these definitions have since then appeared in the literature 
several times. Here are the two expressions, in a slightly different notation from that of the 
original references: 

A = h J dk ^aj c+ ak+ — <4i- a k-j i (3) 

where the Fock space operators a-k^j create and annihilate photons of definite mo- 

mentum k and helicity ±. And, 

A = l - J dr (A(r) ■ H(r) - C(r) • E(r)) , (4) 

where ^E,H^ are the electric and magnetic field operators and ^C,A^ the electric and 
magnetic potential operators [8J. 

All of the expressions d2J), ^ and Q generate the same fundamental symmetry trans- 
formation, albeit in different representation spaces. 

We start the derivation by setting convenient units of eo = fio — 1 for the vacuum electric 
and magnetic constants (thus c = 1 and k = tu). We can then use ^ to write the free space 
Maxwell equations as: 



V x E = iwH H = — i AE 

(5) 

V x H = -iwE E = i AH. 

Equations ^ already reveal that A is an operator that transforms electric fields into mag- 
netic fields and vice versa. Note that invariance under generalized duality transformation 
can be interpreted as equivalence between electric and magnetic responses. In the same way 
that angular momentum generates rotation matrices [H] , let us use A as the generator of a 
continuous transformation parametrized by the real number 6: D{9) = exp(i9A). To obtain 
an explicit expression for the transformation that D(8) performs on the fields, we start by 
showing that A 2 is the identity operator for Maxwell fields. 

A 2 E = A f^j = E, A 2 H = A^j= H, (6) 
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where the equalities in each equation follow from ([5]). Since ^ is valid for all E and H, we 
conclude that A 2 = I for Maxwell fields. Using that A 2 = /, and the Taylor expansion of 
the exponential, the continuous transformation generated by helicity can be written: 

D(0) = exp(i0A) = cos(0) J + % sin(0)A. (7) 

The application of D(6) to both electric and magnetic fields reads 

E g = (cos(0)/ + isin(0)A) E, 

(8) 

H e = (cos(0)/ + isin(0)A) H, 

which, after using ^ again, becomes the well know [151 chap. 6.11] generalized duality 
transformation of electromagnetic fields written in ([!]). 

III. GENERALIZED DUALITY SYMMETRY IN PIECEWISE HOMOGENEOUS 
AND ISOTROPIC MEDIA 

We will now show that generalized duality symmetry can be restored in the macroscopic 
Maxwell's equations independently of the shapes of the material domains involved. The 
macroscopic Maxwell equations are valid whenever the electric and magnetic fields are av- 
eraged over many of the atoms or molecules composing the materials. In this way, for 
most situations, the electromagnetic properties of the materials are determined only by the 
electric permittivities e and magnetic permeabilities \x. 

We consider an inhomogeneous medium Q composed of several material domains with 
arbitrary geometry. We assume that each domain % is homogeneous and isotropic, and fully 
characterized by its electric e, and magnetic constants (we again use eo = f-io = !)• In 
each domain, the constitutive relations are B = /ijH, D = e^E, and the curl equations for 
monochromatic fields read 

V x H = — iwD = — iueiE, V x E = iuB = iw/^H. (9) 

Using A = (l//c)Vx from and to = ko = k/ \feijli we obtain 

AH = -i. /5e, AE = i\ f^H. (10) 

V A*i V e i 

Note that to arrive at this result, the fact that the wavenumber in each medium is k = 
koy/iijli has to be used in the expression of the helicity operator. With this change, we are 
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able to obtain the formal expression of the helicity operator for a material medium, which 
we could not find in the literature. Now, we can normalize the electric field E — > . — E, to 
show that we can recover the exact form of Maxwell's equations in free space Q. Clearly, 
the normalization can only be done when all the different materials have the same ratio 
— = a V i. When the normalization is possible, the electromagnetic field equations on the 
whole medium f2 are invariant under the generalized duality transformations of ([!]). 

The remaining question is what happens at the interfaces between the different domains, 
where the material constants are discontinuous. We now examine the boundary conditions 
in Q. At the interfaces between media, the electromagnetic boundary conditions impose the 
following restrictions on the fields: 



n x (Ei - E 2 ) = 0, n x (Hi - H 2 ) = K, 
n ■ (Di - D 2 ) = ct, n-(Bi-B 2 ) = 0, 



(11) 



where K is the surface current density, a the charge density and n the unit vector perpen- 
dicular to the interface. The boundary conditions can be seen as applying point to point 
to a differential surface area at the interface between the two media [161 chap. 2.8]. Let us 
choose a particular point r on the interface. Assuming no free charges, i.e K = and a = 0, 



equations (11) may be interpreted as a linear transformation applied to the fields at one 



medium which results in the fields at the other medium. Using (11) and the constitutive 
relations, the transformation equation reads: 
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(12) 



z. 



where we have oriented our reference axis so that n 

On the other hand, the generalized duality transformation ([T]) may also be written in 
matrix form: 



Eg 




lcos(0) 


-lsin(0) 
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H, 




lsin(0) 


lcos(0) 




H 



U{9) 
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H 



(13) 
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where / is the 3x3 identity matrix. It is a trivial exercise to check that the transformation 



matrix of (12) commutes with U{6) if and only if ei//^i = €2/^2- in such case, the fields 
in each of the two media can be transformed as in (JT]) while still meeting the boundary 
conditions at point r. We can now vary r to cover all the points of the interface and repeat 
the same argument: The fact that U{9) does not depend on the spatial coordinates allows 
to reorient the reference clXIS clS needed to follow the shape of the interface between two 
media. The derivation is hence independent of the shape of the interface, and we may 
say that the boundary conditions are invariant under generalized duality transformations 
when ex/ fi\ = £2/^2- The above derivations show that both the equations and the boundary 
conditions in f2 are invariant under (jlj) when ^ = aV domain i. As a conclusion, we can state 
that independently of the shapes of each domain, a piecewise homogeneous and isotropic 
system has an electromagnetic response that is invariant under duality transformations if 
and only if all the materials have the same ratio of electric and magnetic constants: 

— = a V domain i. (14) 
Mi 

In this case, since helicity is the generator of generalized duality transformations, the system 
preserves the helicity of light interacting with it. 

Our results are in agreement with Bialynicki-Birula's wave equation for photons propa- 
gating in a linear, time- in dependent, isotropic and inhomogeneous medium. In [T7J §2], he 
shows that the two helicities of the photon are only coupled through the gradient of y 

In the same review [TTl §11], the author discusses the conservation of helicity in arbi- 
trarily curved spacetime. This notable fact is related to the equivalence of the free space 
Maxwell's equations on an arbitrary spacetime geometry and the macroscopic equations on 
a flat spacetime occupied by an anisotropic medium. This equivalence is fundamental in 
transformation optics [18j §4], the theoretical basis for metamaterials. Not all anisotropic 
media represent spacetime geometries, there is a necessary and sufficient condition for it 
[T8l §4]: The electric permittivity and magnetic permeability tensors must be equal to each 
other and equal to a certain function of the spacetime metric = /i lJ = f{g 1 ^) [TSJ §4]. 
At this point, the possible relationship with our results is apparent through the condition 
e lJ = and the conservation of helicity discovered by Bialynicki-Birula. Extending the 
proof in this section to interfaces between anisotropic media could lead to new insights. 



Relation (14) is often referred to as surface impedance matching condition. It has al- 



ready been explored in scattering from spheres [19] (see appendix |Aj) and in the context of 
plasmonic metamaterials [20] ■ The connection with helicity preservation was not considered 
in these references, and, to the best of our knowledge, it has not been considered in any of 



the references where a relationship as the one in (14) is used. The above derivation shows 
that these particular cases are part of a more general, geometry independent symmetry: 
generalized electromagnetic duality. 

In order to illustrate the independence of helicity conservation from geometry we per- 
formed numerical simulations. We analyzed the helicity change (Fig. [I]) for two different 
dielectric structures in free space: A circular cylinder, which is symmetric under rotations 
along its axis, and a curved panfTute like structure without any rotational or translational 
symmetry. Two versions of each structure were simulated, corresponding to two different 
materials: the first one would represent silica by setting e = e g i ass = 2.25 and \i = /i g i aS s = 1- 
In the second material we enforce duality (14), by setting e = fi = e g i ass = 2.25. The inci- 



dent field is a circularly polarized plane wave (i.e. it has well defined helicity) propagating 
parallel to the axis of the cylinder and the curved surface of the panflute. Fig. [T] shows that 



helicity is conserved independently of the spatial symmetries, whenever eq. ( 14 ) is fulfilled, 
i.e. under conditions of duality symmetry. On the other hand, conservation of angular 
momentum, resulting in cylindrical symmetry of the scattered fields, is only achieved in the 
case of the cylinder. 

The most important result of this article is the geometry independent restoration of gen- 
eralized duality symmetry in the macroscopic Maxwell's equations, turning the relationship 
between helicity and duality into a practical tool. Helicity and duality can now be added to 
the toolbox used for the study of light-matter interaction problems within the framework of 
symmetries and conserved quantities. This tool and framework have already been used to 
show that the mechanism of optical spin to orbital angular momentum conversion is an in- 
consistent explanation [21 J. The concept of spin to orbit conversion applied to focusing and 
scattering is masking two completely different physical phenomena related to the breaking 
of two different fundamental symmetries: transverse translational symmetry in focusing and 
electromagnetic duality symmetry in scattering. Appendix [A] contains another example of 
the application of our ideas. The unusual scattering effects from magnetic spheres reported 
in [19] are explained in a straightforward manner with the use of duality, helicity and other 
symmetries and conserved quantities of the system. 



FIG. 1. (Color) Impact of the different symmetries on the field scattered by two dielectric struc- 
tures. The upper row shows the intensity calculations for a symmetric cylinder, while the scatterer 
used for the lower row has a panflute like shape, which lacks all rotational or translational sym- 
metry. In (a) and (c) the structures are made of silica, while in (b) and (d) we enforced duality 
symmetry by setting e = \i = e g i ass = 2.25. The left half side of each subfigure corresponds to the 
scattered field with helicity equal to the incident plane wave A+, while the right half is for the 
opposite helicity A_. In order to use the same color scale for all figures, the right half side is mul- 
tiplied by the factor in the upper right corner. As can be seen, the (lack of) cylindrical symmetry 
of the structures results in (non-)cylindrically symmetric field patterns, which is consistent with 
the geometry of each case. On the other hand, both scatterers behave identically with respect to 
conservation of helicity, which is seen to depend exclusively on the electromagnetic properties of the 
material. In further simulations we verified that helicity conservation is independent of the angle 
of incidence. This figure is a clear illustration that cylindrical symmetry is related to conservation 
of angular momentum, while duality is related to conservation of helicity. 

IV. PREPARATION, MEASUREMENT AND ANALYSIS OF MAXWELL FIELDS 
WITH WELL DEFINED HELICITY 

Let us now turn to the practical matter of how to prepare and measure electromag- 
netic helicity in optical experiments. To the best of our knowledge, there has not been any 
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proposal for controlling and analyzing the helicity content of electromagnetic beams. As 
discussed in the Introduction, the helicity operator provides information about the polar- 
ization of the plane wave modes composing the electromagnetic field. From a purely optical 
perspective, it seems intuitive to use a combination of lenses, waveplates and polarizers to 
study helicity. This is indeed our approach as can be seen in Fig. |2^c). The detailed analysis 
is somewhat involved since it must be done for general (non-paraxial) Maxwell modes. For 
this reason, we provide the theoretical background needed to understand the structure of 
a useful set of modes of well defined helicity, and we explain the use of aplanatic lenses 
and waveplates to control and measure helicity. Additionally, we identify the experimental 
signature of a scatterer that preserves angular momentum but partially converts helicity. 
All this background will allow the analysis of our experimental setup and results using sym- 
metries and conserved (or non-conserved) quantities. In our experiments, we illuminate 
cylindrically symmetric nanoapertures with beams of well defined helicity and analyze the 
forward scattered field in terms of its helicity content. 

A. Electromagnetic modes with well defined angular momentum (J z ) and helicity 
(A): Bessel beams. 

Typically, modes of well defined helicity, when studied in real space coordinates, contain 
all the possible polarizations. Let us take a look at a set modes which are very relevant in 
our experimental setup: Bessel beams. Bessel beams are cylindrically symmetric transverse 
(i.e. zero divergence) solutions of Maxwell's equations. They have well defined values of 
energy H and third components of angular and linear momentum J Z ,P Z . Among the several 
families of Bessel beams, we will use the one in which they are also eigenstates of helicity, 
A. Helicity commutes with all the other operators. The set of observables H, J Z ,P Z and A 
fully determines the electromagnetic field. This kind of Bessel beams have been reported, 
for example in [22J. A complete derivation can be found in [211 app. B]. 

Below we give the expressions for the two types of vector wave functions C mpz and D mpz 
having sharp values of H = k, P z = p z , J z = m and A equal to —1 and +1 respectively. We 
use cylindrical coordinates [p, 9, z] for the spatial variables and the helical basis [f , 1, z] for 
the vectorial character of the fields, where 1 = (x + iy) / \/2, f = (x — iy) / \/2. An implicit 
harmonic exp(— iwt) dependence is assumed. 
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C mpz (p, 0, z) = A(z) exp(im9) B + J m+1 (p p p) exp(id)r + £_J m _i(p p p) exp(—iO)l + iV2p p J m (p p p)z 
D mpz ((3, 9, z) = A(z) exp(im6>) B_J m+1 (p p p) exp(i<9)r + B + J m -i{p P p) exp(-i0)l - iV2p p J m (p p p)z 

(15) 

where p 2 p = k 2 — p 2 z = pi + p 2 ,, J m {-) are the Bessel functions of the first kind, the am- 
plitude A(z) = \J^i m e x p(Wzz)^/k, and B± = (k ±p z ). These modes form a complete 
orthonormal basis of transverse Maxwell fields. 

Note that the structure of the modes with opposite helicity is very similar and it is only 
distinguished by constants multiplying the different vector components. In the collimated 
limit, when ^ — * (Pz ^ k so that B + — > 2 and £>_ — > 0), both C mpz and D mpz ap- 
proach pure right circularly (RC) and left circularly (LC) polarized modes respectively. The 
other polarization components, the opposite circular and the longitudinal z component, are 



strongly attenuated in the collimated regime. Then, as — > 0, the modes in (15) can be 
approximately written as: 

C mPz (p,e,z) w x r^i m+1 e xp(t(p z z))J m+1 (p p p)ex V (t6(m + l))r, (16) 

V mpz (p,6,z) w x /^i m+1 exp(t(p z z))J m ^(p p p)exp(z6(m - 1))1. (17) 



7T 

The error in the intensity that we are making with this approximation is of the order (p p /k) 2 . 
This asymptotic property will allow us to prepare and analyze light beams with well defined 
helicity using simple optical elements. 

B. Preparation and measurement of Maxwell fields with well defined helicity 



Relations (16) and (17) indicate that for collimated light beams, controlling the polar- 
ization of the field allows us to control its helicity content. For example, let us consider our 
experimental setup of Fig. [2|c). We start with a collimated Gaussian beam with diameter 
w = 5 mm. We use a linear polarizer and a quarter waveplate after the laser source in 
order to obtain a collimated LC polarized Gaussian beam. If we expand such beam with 



the collimated modes of Eqs. (16) and (17), the D mpz components would dominate because 
of the polarization selectivity of the linear polarizer followed by the quarter waveplate. The 
LC/RC polarization intensity ratio of a Gaussian beam with well defined A = +1 will be of 
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the order (wk)~ A , which in our case would be of the order of 10~ 19 . A comparison of this 
figure with the typical extinction ratios of commercial polarizers, 10 -5 , indicates that the 
polarizer is the limiting factor when preparing a collimated beam with well defined helicity. 

After the polarizer and waveplate, we have prepared an electromagnetic field which is 
a superposition of only D-type Bessel beams. The amplitudes of the D mpz for different m 
and p z will be given by the shape of our beam. In particular, if the collimated beam is 
cylindrically symmetric, the value of m controls the azimuthal phase of the field, as seen in 



(17). Since our LC Gaussian beam does not contain any azimuthally varying phase, only 
terms with m — 1 are possible in its expansion. We have hence prepared a J z — +1, A = +1 
light beam. Similarly, in order to project a collimated beam onto states of well defined 
A = ±1 it is possible to use a quarter waveplate and a linear polarizer in two orthogonal 
settings, as we do in our experiments. 

While preparation and measurement are performed on collimated beams, where the ap- 



proximations (16) and (17) hold, the actual interaction of light with the target may happen 
in between two microscope objectives. For example, the first one focuses the incident beam 
onto the target from one side and the second one collects and collimates the output scattered 
field from the other side. Assuming that the two microscope objectives work as aplanatic 
lenses (23], they will not affect the helicity state of the light beam [21], and, if perfectly 
aligned with the incident beam will not affect its angular momentum J z either: a formal 
proof of these assertions can be found in [2H app. C]. Since a lens does not change the 
energy either, the lenses in our experimental set-up only redistribute the amplitudes of the 
different p z modes. This has two important consequences. First, by using simple prepara- 
tion methods it is possible to illuminate a target with a strongly focused beam possessing 
well defined A, and second, that whatever helicity change is observed by the measurement 
apparatus can be attributed to the interaction of the light beam with the target. 



C. Conservation of angular momentum and helicity change 

Let us assume that the target in question is cylindrically symmetric, that is, it preserves 
J z . Whether it preserves helicity or not depends only on the material properties, as illus- 
trated in Fig. [Tj The transformation of a light beam interacting with such a target can be 
represented by a transfer matrix between Bessel modes (C mp2 , D mp J. Since the transfor- 
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mation leaves the angular momentum invariant, the transfer between Bessel modes is only 
allowed between modes with the same index m. If the transformation leaves the helicity 
invariant, as in the case of the aplanatic lens, C-type modes will be transferred to C-type 
modes, and similarly for D-type modes. If a cylindrically symmetric target partially con- 
verts the helicity of an incident beam, there will be a transfer between C and D modes with 
the same index m. Let us study such case in detail. 

Let us say that we start with a J z — +1, A = +1 beam which we focus on a cylindrically 
symmetric target. If a portion of this beam undergoes an angular momentum conserving 
helicity transformation, becoming J z — +1, A = — 1, then after collimation, the portion of 
the beam which has undergone the change can be expanded as a sum of only modes of the 
type 

/ D 

C lpz (p,d,z) « J— i 2 exp(i(p z z))J 2 (p p p) exp(i29))i, (18) 

while the portion which did not experience a helicity change will necessarily be a sum of 
only modes of the type: 

IT) a 

Di/>,(P,M) ~ \ — i 2 exp(i(p z z))J (p p p))\. (19) 
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Expressions (18) and (19) have been obtained by setting m — 1 in (16) and (17). We see 
that the angular momentum preserving helicity transfer always leaves an azimuthal phase 
imprint in the collimated regime. In our case, an original LC polarized beam with no phase 
singularity transforms into a RC polarized beam with a phase singularity of order two and 
its corresponding zero of intensity at p = ( Jl8| ) (Note that Ji(0) = V I ^ 0). This vortex of 
charge two is hence a signature of an angular momentum preserving and helicity changing 
scattering. A projective measurement onto the helicity value opposite to the input one 
followed by a Charge Couple Device (CCD) camera can be used to look for such a signature. 



V. HELICITY TRANSFORMATIONS IN NANOAPERTURES 

In order to investigate helicity changes in light matter interactions, we performed a series 
of experiments with nanoapertures in a gold film over a glass substrate. The transmission 
of light through isolated nanoholes was first studied in the seminal paper of Bethe [25] 
and is now understood [26] to be crucially affected by the effect of localized plasmons and 
surface modes at the metal-dielectric interfaces. Applications of this kind of nanostructures 
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include optical trapping [27], funneling of light [28] and reshaping of optical fields |29j. In 
our experiments we observed that helicity is transformed, even though the system we probe 
is cylindrically symmetric. We thus experimentally verify that helicity is decoupled from 
angular momentum. We analyze the complete set-up and measurement results only in terms 
of symmetries and conserved quantities. This methodology allows us to make qualitative and 
quantitative predictions and, as a result, identify the exact reason for the relatively large 
helicity conversion in the nanoapertures: The coupling of light to surface modes, which, 
being linear combinations of modes of opposite helicity, strongly break duality symmetry. 

Figure [2]^a) sketches the interaction of light and a scatterer consisting of a circular 
nanoaperture milled in a gold film. Two symmetries of this simple scatterer are notice- 
able: rotational symmetry around the propagation direction, i.e. the z-axis, and mirror 
symmetry upon reflection through any plane containing the z-axis. Rotational symmetry 
along the z-axis implies conservation of the component of angular momentum projected 
along such axis, i.e. J z = z • J. All the machinery for the preparation, measurement and 



analysis of helicity, explicitly developed in section |IV A| for cylindrically symmetric targets, 
is fully applicable to our experimental setup. 



A. Experimental setup, methods and measurements 

The system we have experimentally tested consists of a set of isolated nanoapertures of 
different sizes, which were milled with a focused ion beam in a gold layer of 200 nm, deposited 
on top of a 1 mm thick glass substrate. The diameters of the nanoapertures ranged from 
150 to 580 nm. Using a Scanning Electron Microscope (SEM), we carefully measured the 
diameters of all the apertures and checked that all the apertures had aspect ratios between 
0.95 and 1. We probed the nanoapertures with a continuous wave laser with a wavelength of 
Ao = 633 nm. The preparation of the probing beam was done as follows. First, we collimated 
the laser beam and then we used a set of linear polarizers and waveplates to ensure a left 



circularly polarized light beam. As explained in section IV B, when this collimated field is 
decomposed in modes of well defined J z and A, the components with J z = 1, A = 1 are 
overwhelmingly dominant. This collimated helicity field was subsequently focused with a 
microscope objective with a numerical aperture of NA=0.5. Since the transformation of an 
aplanatic lens preserves helicity, we were able to generate a focused electromagnetic field 
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FIG. 2. (Color) (a), Sketch of the symmetries of the system we are probing. Light with a well 
dehned helicity, represented by a red propeller, impinges on a cylindrically symmetric aperture in 
a thin metallic film. The output light is analyzed in terms of its helicity content, represented as 
red propellers when it is the same as the incident helicity, or blue ones with opposite handedness 
when the helicity is the opposite to the incident one. Note that helicity always flips upon a mirror 
transformation. This is graphically shown in the sketch by the addition of a mirror plane, (b), 
Scanning Electron Microscope (SEM) image of the apertures that we tested with our system. We 
probed a total of 212 different isolated apertures. The separation between apertures was either 5 
lim or 50 //m and we could not find relevant differences between the two sets, (c), Experimental 
set-up. An incoming collimated beam is circularly polarized with a set of waveplates and focused 
to address the isolated nanohole, which is centered with respect to the beam with a nanopositioner. 
The transmitted light is then collected and analyzed with another set of waveplates and a Charge 
Coupled Device (CCD) camera. 

with a well defined helicity. The focused field was then allowed to interact with one of the 
isolated nanoapertures. We carefully positioned the nanoaperture on the symmetry axis of 
our optical system by means of a set of piezo-stages. Subsequently, the scattered light was 
collected and collimated with another microscope objective of NA=0.9. Once again, this 
lens did not affect the helicity of the beam, and as such, after collimation, we were able 
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FIG. 3. (Color) Helicity transformations through nanoapertures. (a), Dependence of the ratio of 
transmitted helicities on the aperture size. All points correspond to highly symmetrical apertures, 
the sizes of which were measured with an SEM. (b), Numerically calculated spatial pattern for the 
direct and transformed helicity fields transmitted through a cylindrical aperture of 300 nm 

to analyze the helicity with another set of waveplates and polarizers, obtaining two very 
different spatial profiles for fields with J z — 1, A = 1 and J z — 1, A = — 1. The light was 
detected with a CCD camera. 

Our results show that there is always a helicity transformation in the transmitted light. 
This can be seen in Fig. [3|a), where we plot the power ratio between the two transmitted 
helicities, 7, as a function of the aperture sizes. The smallest conversion we measured was for 
the largest holes, 7580 = 0.08 ±0.02. In contrast, the helicity transformation measured in the 
same sample through the glass alone, corresponding to an infinite aperture, is 7; n f ~ 10~ 3 . 
This conversion value can be due to left over gold nanoparticles distributed over the glass 
surface. Also, the helicity transformation by the focusing and collimating lenses alone was 
even smaller, of the order of 7i ens ~ 10 -4 , which is consistent with the fact that perfect 
aplanatic lenses should preserve helicity. In Fig. |3^b), we display the typical spatial patterns 
for the two output helicities of the light scattered from a perfect cylindrical aperture, as 
calculated with a semi-analytical method [30]. We numerically checked that this output 
conserves angular momentum but, as can be seen, breaks helicity conservation (duality 
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FIG. 4. (Color) Projective measurement of helicity. Upper (lower) row shows the results for an 
aperture of 500 nm (300 nm) . (a) and (d) , Experimental results showing the transmitted light with 
helicity identical to the incident light, (b) and (e), Experimental results of transmitted light with 
opposite helicity. The cylindrical asymmetries are due to the finite rejection power of the polarizer. 
This is clearly seen with the numerically simulated patterns, (e) and (f), where the asymmetry 
appears only after including the experimental parameters of the used waveplates. 

symmetry). The spatial shape and optical vortex content of beams that have undergone an 



angular momentum preserving but helicity flipping interaction is explained in section IV C 



In order to test that our experimental results are consistent with helicity changes with 
conservation of angular momentum, we analyze the CCD images. In Fig. [4] we show typi- 
cal experimental results and their comparison with numerical calculations for two different 
aperture sizes. In the left column (Figs. Qa) and (d)) we show the components of the out- 
put field with the same helicity as the input, A + . The observed field pattern is a typical 
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Airy pattern arising from the subwavelength dimensions of the nanoaperture and the finite 
numerical aperture of the collection microscope objective, as expected from Fig. ^b). On 
the other hand, the central column (Fig. [4](b) and (e)) shows the field with opposite helicity, 
A_. The absence of singularities in Figs. |4|a) and (d), and the presence of two singularities 



in Fig. Kb) and (e), are fully consistent with the results of section IV C 



The differences between our experimental results and the ideal case of Fig. |3jb) for the 
helicity transformed transmission are due to the finite extinction ratios of our polarizers. 
We now discuss this assertion. First of all, electromagnetic modes with well defined J z 



have to be cylindrically symmetric. Additionally, as we show in section IV B, collimated 
beams with well defined angular momentum and helicity always present a phase singularity 
of order q = J z — A in the dominant polarization, i.e. the wavefront of the field is twisted 
around the center q times 2tc. In our case, this translates to a smooth phase front, with no 
phase singularity, for the directly transmitted helicity (q = 1 — 1 = 0) and a second order 
singularity for the helicity transformed mode (q = 1 — (— 1) = 2); which is consistent with 
the experimental images of Fig|4} In practice, we could not avoid a small leakage from the 
direct helicity to the transformed helicity. As a result of such a superposition, the second 
order singularity splits into two singularities of order 1. Thus, the intensity pattern is no 
longer cylindrically symmetric. In order to prove this point, we show in Figs. |4]^c) and (f) 
the coherent superposition of numerically calculated images of unmixed modes (as those of 
Fig. |3^b)), with relative amplitudes given by the extinction ratios of our set of polarizers, 
5 x 10~ 5 . Given that the leakage through the waveplates is the same for all apertures, one 
would expect that its effect should be smaller for those scatterers with a higher 7. This is 
indeed the case seen in Fig. |4j which shows a larger spatial separation of the singularities for 
the larger aperture. The helicity transfer for the 500 nm aperture, 7500 = 0.07 ± 0.01, is a 
factor 2.3 smaller than for the 300 nm aperture, 7300 = 0.16 ± 0.03. The only remaining free 
parameter is the relative phase between the two modes, the effect of which is to rotate the 
whole pattern around its center. By comparing the experimental results with the calculated 
pattern we can infer that the relative phase between the modes was tt/3, for the 500 nm 
aperture (Fig. [l|c)), and for the 300 nm (Fig. |i](f)). We conclude that our measurements 
are consistent with the fact that in our system, the angular momentum is conserved, but 
helicity is not. 
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B. Analysis of the experimental measurements 



According to the ideas presented in this paper, the observed helicity change implies that 
electromagnetic duality is broken in our samples. In order to identify the exact mechanism 
of duality breaking, it is important to first consider the multilayer system air-glass-gold-air 



without the nanoaperture. Considering condition (14), duality is obviously broken by just 
the multilayer alone, but the helicity transfer in the absence of the nanoaperture has been 
numerically shown to be around 10~ 6 for our collection objective. This value is one order 
of magnitude smaller than the typical measurement noise, which is related to the polarizer 
extinction value of 5 x 10~ 5 . The experimental observation of much higher transformation 
ratios must hence be tied to the nanoapertures. 

The transmission and reflection in a multilayer system is best studied using plane waves. 
For a single plane wave with momentum k, the two helicity states are the two states of 
circular polarization, and can be obtained by linear combination of its s (transverse electric) 
and p (transverse magnetic) components: sip. See (211 a PP- A] for a general derivation 
of this relationship, which also applies to multipolar fields and Bessel beams. Different 
scattering coefficients for s and p will hence mix the two helicity modes. This idea has been 
recently applied to the analysis of resonances in spheres [31] . If the particular multilayer 
presents any resonance, for either s or p, the helicity transfer will be enhanced in its vicinity: 
A pure s or pure p mode is an equal weight combination of the two helicities, and hence 
strongly breaks helicity conservation (duality symmetry). Our system indeed presents two 
resonances for non-propagating modes. These two resonances are related to surface modes on 
the metallic surfaces. It is well known that through the scattering of the nanoaperture, the 
incident field can excite surface plasmon polaritions (SPP) at the interface. Since SPPs are 
p-polarized waves [32], this produces an asymmetric response of the SPP with regard to the 
transmitted s and p-polarized components. This additional SPP induced electromagnetic 
asymmetry dramatically enhances the helicity transfer, even in the propagating modes, 
and makes it experimentally detectable. According to this explanation, the helicity transfer 
should increase for modes in the proximity of the plasmon resonance, i.e. for large transversal 
momenta. This actually explains the trend in Fig. [2] that smaller holes present a larger 7 
value: smaller holes have a higher coupling to large transversal momenta, and in particular to 
the surface modes. We can then conclude that the most important role of the nanoaperture 
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in the helicity transfer is to couple the incident field to the resonances and, in particular, to 
the surface plasmon polaritons. 

Very similar experimental setups and measurements have been previously reported. Start- 
ing with the work of Gorodetski et al. [33J, we can also find [M] and very recently [35]. In all 
these cases, the results have been analyzed by means of the spin (S) to orbital (L) angular 
momentum conversion mechanism. As already mentioned, [21] shows the inconsistency of 
the spin to orbital explanation, whose root is the fact that only the components of J = S + L 
are the generators of rotations for the electromagnetic field. The components of L and S, 
when separately considered, do not generate meaningful symmetry transformations of the 
electromagnetic field. This is reflected by the fact that the application of the operator cor- 
responding to any of their components, for instance L z or S z , breaks the transversality of 
the fields [20], [33 page 50]. In this paper, we have used the generator of rotations along 
the z axis J z , and the generator of generalized duality transformations A, to analyze our 
experimental results from the point of view of symmetries and conserved quantities. This 
framework has allowed us to identify the root cause of the relatively high helicity change 
by means of qualitative and quantitative considerations. Spin to orbit angular momen- 
tum conversion, due to its inconsistent nature, does not allow this level of understanding. 
The analysis contained in this section explains the experimental observations of the above 
given references, including some of the cases in [33] and [34J where the apertures were non- 
cylindrical. 



VI. CONCLUSION 

In this article, we have shown that the restoration of generalized duality symmetry is 
possible for the macroscopic Maxwell's equations, even though the microscopic equations 
are rendered asymmetric by the empirical absence of magnetic charges. The restoration of 
the symmetry is independent of the geometry of the problem: a system made of piecewise 
isotropic and homogeneous domains of different materials characterized by electric and mag- 
netic constants (ei,/ij) is invariant under generalized duality transformations if and only if 
— = a V media i. This result is independent of the shapes of the domains. With this result, 
the known relationship between helicity and generalized duality transformations, namely 
that the former is the generator of the latter, is turned into a simple and powerful tool for 
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the practical study of light-matter interactions using symmetries and conserved quantities. 
Armed with it, we have experimentally investigated helicity transformations in focused light 
fields that interact with cylindrical nanoapertures in a gold film over a glass substrate. An- 
alyzing the results by means of symmetries and conserved quantities, including duality and 
helicity, we have been able to conclude that the role of the nanoapertures is to allow light to 
couple to resonant modes of the system where duality is strongly broken: this is what renders 
the helicity transformation effect observable in our experimental set-up. In our experiments, 
we have used a fact of crucial practical importance: Measurement and preparation of beams 
with well defined helicity can be done with very simple optical elements. 

We are confident that we are proposing a useful tool. In the article, we have shown 
how to apply it to an experiment and gain valuable insight with it. In appendix [A] we use 
it to explain some unsual scattering effects in a straightforward way. In [21J, it allowed 
to proof the inconsistency of the concept of optical spin to orbital angular momentum 
conversion in focusing and scattering, and to propose a substitute framework based on 
helicity. Additionally, we are preparing a manuscript were the role of duality symmetry in 
molecular optical activity is identified. 

Our results may be useful in other fields. For example, they may prove important in the 
field of metamaterials and transformation optics |38j, which is dramatically extending the 
range of wavelengths where effective electric and magnetic constants can be engineered. The 
transfer of helicity between light and matter remains an open line of research, which could 
have importance in the fields of plasmonics and "spintronics" [39], where the control of the 
helicity of electrons is crucial. Finally, it can be seen that the same tools we have developed 
here can be successfully used to explain recently reported effects in electron beams jlO]. This 
parallelism is an encouraging sign towards the possibility of simulating particle interactions 
on an optical table pETj . 

Appendix A: Scattering effects for magnetic spheres 

Consider the unusual scattering effects for magnetic spheres reported by Kerker |19j . One 
of them refers to the fact that a plane wave impinging on a vacuum embedded sphere with 
- = 1 does not produce any backscattered field (at a 180 degrees scattering angle). This 
effect, which has been referred to as an anomaly [12], can be easily understood using our 
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results. Let us take as incident field a circularly polarized plane wave with is momentum 
aligned along the z axis. Its angular momentum is also aligned with the z axis and, in 
natural units of h = 1, equal to ±1 depending on the handedness of the polarization. We now 



know (equation (14)) that, based only on the properties of the materials, helicity has to be 
preserved in the interaction between the plane wave and a dual sphere. Let us now assume 
that there exist a component of backscattered field at 180 degrees, that is, a plane wave 
whose linear momentum is the negative of the linear momentum of the incident plane wave. 
Recalling that A = J • P/|P|, we see that, to preserve helicity, the angular momentum of 
the backscattered plane wave must also change sign with respect to the angular momentum 
of the incident plane wave. But such change is impossible: the rotational symmetry of 
the sphere implies that angular momentum is preserved in all axes, in particular along 
the axes shared by the incident and the backscattered plane waves. As a consequence, 
the backscattering amplitude must be zero. Since this argument applies independently to 
both circular polarizations, the backscattering gain will be zero for any polarization of the 
incoming plane wave. What is sometimes referred to as the first Kerker condition is hereby 
explained. Interestingly, zero backscattering from dual objects has already been the object 
of investigation by the radar community jl3] : in these works, the connection between duality 
and helicity is not recognized. 



In the same paper, Kerker finds that upon scattering off a vacuum embedded sphere 
with - = 1, the state of polarization of light is preserved independently of the scattering 
angle. The root cause of such interesting phenomenon is the simultaneous invariance of the 
system with respect to generalized duality transformations, due to the materials, and any 
mirror operations through planes containing the origin of coordinates, due to the geometry. 
In the helicity basis, the 2x2 scattering matrix between an incident and a scattered plane 
wave jH chap. 3] must be diagonal because of helicity preservation. Additionally, it must 
also preserve the linear polarizations parallel and perpendicular to the plane containing the 
two plane wave momentum vectors, because a mirror operation across such plane leaves the 
sphere and both momentum vectors invariant. Using then that helicity flips with mirror 
operations, it can be easily shown that all the 2x2 scattering matrices are indeed diagonal 
and hence preserve the state of polarization between any pair of incident and scattered plane 
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waves. 
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